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Size  Dependent  Meehanieal  Behavior  of  Free-Standing  Glassy  Polymer  Thin  Films 

ABSTRACT 

Meehanieal  properties  of  nanoseale  free-standing  polymer  thin  films  exhibit  size-dependenee  due  to  surfaee  effeets. 
However,  it  has  remained  a  ehallenge  to  assoeiate  the  length  scales  where  such  differences  emerge  with  bulk  polymer 
properties.  Here  we  utilize  molecular  dynamics  simulations  to  uncover  the  dependence  of  elastic  modulus  of  free¬ 
standing  films  on  film  thickness  and  bulk  properties.  Comparison  of  glass  transition  temperature  (Tg)  and  modulus 
(E)  indicates  that  Tg  converges  to  bulk  value  slightly  faster  as  film  thickness  increases.  The  free  surface  effects  that 
give  rise  to  a  depression  in  E  and  Tg  are  observed  to  be  stronger  for  polymers  with  weaker  intermolecular 
interactions.  Our  simulations  suggest  that  the  length  scale  of  perturbation  of  free  surface  is  only  about  several 
nanometers,  but  the  effect  is  large  enough  that  only  films  of  100  nm  or  larger  exhibit  negligible  surface  effects. 
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ABSTRACT 

Mechanical  properties  of  nanoscale  free-standing  polymer  thin  films  exhibit  size-dependence  due 
to  surface  effects.  However,  it  has  remained  a  challenge  to  associate  the  length  scales  where  such 
differences  emerge  with  bulk  polymer  properties.  Here  we  utilize  molecular  dynamics 
simulations  to  uncover  the  dependence  of  elastic  modulus  (E)  of  free-standing  films  on  film 
thickness  and  bulk  properties.  Comparison  of  glass  transition  temperature  {Tg)  and  E  indicates 
that  Tg  converges  to  bulk  value  slightly  faster  as  film  thickness  increases.  The  free  surface  effects 
that  give  rise  to  a  depression  in  E  and  Tg  are  observed  to  be  stronger  for  polymers  with  weaker 
intermolecular  interactions.  The  most  intriguing  aspect  of  our  study  is  the  finding  that  despite 
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the  observed  deerease  in  modulus  of  film  up  to  a  thiekness  over  100  nm,  the  loeal  stress 
distribution  reveals  that  the  preserved  length  seale  of  perturbation  of  free  surfaee  is  only  about 
several  nanometers.  {Abstract  Figure:  Fig.  1 ) 

INTRODUCTION 

Ultra-thin  polymer  films  with  thickness  in  the  range  of  tens  of  nanometers  are  becoming 
relevant  for  a  wide  array  of  applications  in  the  field  of  nanotechnology  and  bioengineering,  such 
as  nanoelectronics  \  nanocomposites  coatings  polymeric  membranes  and  biosensors  For 
instance,  in  nanoelectronics,  the  thermal  and  mechanical  responses  of  nanopattemed  polymeric 
structures  obtained  from  lithography  is  important  in  determining  the  stability  and  reliability  of 
these  systems  In  thin  supported  and  free-standing  films,  the  effects  of  free  surfaces  on  the 
molecular  structure  and  dynamics  of  polymer  chains  fundamentally  changes  the  behaviors  of 
these  systems  compared  to  bulk,  and  these  effects  become  significantly  more  pronounced  as  the 
film  thickness  reaches  below  100  nm.  The  free  surface  effects  are  complex  and  depend  on  many 
factors  such  as  molecular  weight  of  polymers,  thickness  of  the  film,  and  polymer  chemistry,  but 
a  common  observation  is  that  the  presence  of  surfaces  in  polymer  thin  films  may  result  in  a 
notable  reduction  of  thermo-mechanical  properties  by  shifting  the  relaxation  time  and  enhancing 
the  molecular  mobility  Understanding  fundamental  mechanisms  pertaining  to  the  thermo¬ 
mechanical  properties,  in  particular,  the  glass  transition  temperature  (Tg)  and  elastic  modulus  (E) 
of  sub- 100  nm  polymer  thin  films  are  important  for  numerous  technological  applications  of  soft 
nanomaterials,  and  are  being  studied  through  novel  experimental  and  computational  techniques 
suited  for  the  length  scales  of  interest. 
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Previous  studies  have  shown  that  intermolecular  noneovalent  weak  interaetions  between 
polymer  ehains  are  erueial  in  determining  thermo-meehanieal  responses  of  amorphous  polymers. 
Experimental  and  theoretieal  studies  have  eonsistently  shown  that  intermoleeular  forees  between 
polymer  ehains  ean  be  tailored  to  tune  Tg  of  bulk  amorphous  polymer  Concerning  the 
mechanical  behaviors,  the  investigations  of  energy  partitioning  have  shown  that  the  elastic 
modulus  and  yield  response  of  amorphous  polymer  are  mainly  dominated  by  intermolecular 
noneovalent  interactions  (i.e.  van  der  Waals  interaction)  and  these  mechanical  properties 
are  intimately  related  to  segmental  dynamics  of  chains  and  physical  aging  phenomena  Low 
molecular  weight  plasticizing  and  antiplasticizing  additives  are  often  used  to  soften  or  stiffen 
polymer  by  changing  their  free  volume  and  molecular  packing  Experimental  investigations 
and  simulations  have  suggested  that  the  interface  and  free  surface  effects  on  thermo-mechanical 
behaviors  of  polymer  films  can  be  influenced  by  numerous  factors  associated  with 
intermolecular  forces,  such  as  retained  solvents  ’  ,  monomeric  interactions  ’  ’  and  chain- 
tacticity-related  density  variations 

More  recent  efforts  have  been  engaged  in  understanding  the  molecular  origins  of  effects  of  free 
surfaces  on  thermo-mechanical  properties  of  polymer  films.  The  thickness-dependent  behaviors 
of  glass  transition  temperature  and  elastic  modulus  of  polymer  films  have  been  characterized  via 
experiments  and  simulations.  Below  Tg,  the  elastic  moduli  of  polystyrene  (PS)  and  poly(methyl 
methacrylate)  (PMMA)  films  and  nanostructures  decrease  dramatically  below  -40  nm  measured 
by  a  buckling-based  metrology  Torres  et  al.  show  that  the  Tg  and  elastic  moduli  of  sub-100 
nm  methacrylate  polymer  films  both  decrease  with  film  thickness,  but  the  decrease  in  modulus  of 
films  is  not  correlated  with  the  observed  decrease  in  7^  A  study  by  Declambre  et  al.  shows 
that  the  elastic  modulus  of  free-standing  antiplasticized  PMMA  nanostructure  decreases  below 
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bulk  value  when  the  width  is  below  50  nm  measured  by  beam  bending  deformation  On  the 
other  hand,  reeent  experimental  evidenees  show  surfaee  stiffening  as  observed  in  polymer  thin 
films  in  the  rubbery  regime  above  Tg  in  whieh  the  underlying  physies  of  meehanieal 
behaviors  eould  be  fundamentally  different  from  that  below  Tg. 

However,  due  to  the  diffieulty  in  measuring  meehanieal  responses  of  polymer  free-standing 
thin  film  with  free  surfaees  only  (no  substrate)  as  illustrated  in  Figure  1  (a),  current  experimental 
techniques  of  measuring  elastic  modulus  of  film  inevitably  introduce  substrate  and  retained 
solvents  effects,  which  have  been  recently  found  to  be  crucial  in  determining  the  surface 
dynamics  at  nanoscale  More  importantly,  how  intermolecular  forces  between  polymer  chains 
interplay  with  free  surface  that  changes  influence  in  Tg  and  mechanical  responses  of  polymer 
films  remains  challenging  to  discern,  and  the  physics  of  this  intriguing  shift  is  yet  to  be  fully 
explored.  Nevertheless,  direct  comparison  of  polymers  with  different  chemistry  via  experiments 
and  all-atomistic  simulations  remains  a  challenge  since  these  parameters  cannot  be  continuously 
and  independently  varied.  Simplified  potentials  in  coarse-grained  molecular  dynamics  (CGMD) 
simulations  of  free-standing  film  with  free  surfaces  provide  a  clear  benefit  in  this  regard  in 
addition  to  computational  efficiency. 

To  facilitate  a  deeper  understanding  of  the  interplay  between  intermolecular  forces  and  free 
surface  effects  on  changing  the  thermo-mechanical  responses  of  polymer  thin  films,  we  employ 
coarse-grained  molecular  dynamics  simulations  to  investigate  independent  effects  of  these 
parameters.  For  this  purpose,  we  establish  a  CGMD  approach  that  builds  on  realistic  parameters 
derived  from  all-atomistic  MD  simulation  of  PMMA.  By  employing  this  CG  model  as  a  starting 
point,  we  carry  out  large-scale  MD  simulations  to  assess  the  free  surface  effects  by  examining 
thickness-dependence  behaviors  of  7^  and  modulus  of  free-standing  thin  films  with  two  free 
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surfaces.  For  modulus  of  polymer,  we  focus  only  on  the  elastic  region  at  a  glassy  state  below  Tg, 
where  polymers  are  deformed  signifieantly  faster  than  the  intrinsie  relaxation  time  and  their 
mechanieal  properties  are  almost  independent  of  temperature.  The  interplay  between 
intermolecular  forces  and  free  surface  effect  on  film  Tg  and  modulus  is  evaluated  by  varying  the 
depth  of  the  effective  weak  interaction  potential  between  monomers.  The  interfacial  and 
characteristic  length  scales  that  eorrelate  with  the  scales  where  Tg  and  E  of  thin  film  deviate  from 
the  bulk  are  examined  and  eompared.  A  bilayer  eomposite  model  derived  from  the  spatial  stress 
distribution  is  used  to  explain  the  physieal  origin  of  the  size  effect  on  film  modulus. 

METHODS 

Coarse-grained  Model  Development 

In  order  to  simulate  films  with  a  characteristie  dimension  (i.e.  thiekness)  reaching  -'100  nm, 
we  first  present  a  coarse  grained  (CG)  modeling  approach  based  on  the  output  of  all-atomistie 
simulations  of  PMMA.  The  one  bead  per  repeating  unit  mapping  from  the  all-atomistic  model  to 
CG  model  and  the  resultant  CG  free-standing  film  configuration  are  illustrated  in  Figure  1  (b). 
The  bonded  interaetion  terms  including  the  bond  stretch  and  angle  bending  terms  of  the  CG 
model  are  derived  to  mateh  the  corresponding  probability  distribution  functions  of  all-atom 
model.  For  the  intermolecular  noncovalent  interactions  between  ehains  (ineluding  interehain  and 
intrachain  interactions),  we  employ  the  truneated  Lennard- Jones  (LJ)  12-6  potential: 

U^(r)  =  4s  -fyj  /•<«,  (1) 

where  e  is  the  depth  of  the  potential  well,  and  a  eorresponds  to  the  radial  distance  where  the 
potential  energy  is  zero.  These  two  parameters  are  ealibrated  to  mateh  the  experimental  density 
at  300  K  and  Tg  of  bulk  PMMA.  The  resultant  values  of  a  and  s  are  6.5  A  and  0.3  kcahmol. 


respectively,  which  yield  a  density  of  1.17  g/cm^  at  300  K  and  a  bulk  Tg  of  388.6  K  for  our  CG 
model.  The  cutoff  distance  of  the  intermolecular  interaction  is  chosen  to  be  15  A.  Detailed 
description  of  the  coarse-graining  strategy  and  simulation  protocols  can  be  found  in  our  previous 
work 

Utilizing  the  atomistically  informed  CG  model  as  the  generic  polymer  model,  we  are  able  to 
explore  the  effect  of  intermolecular  noncovalent  interactions  between  chains  on  the  thermo¬ 
mechanical  properties  (i.e.  Tg  and  elastic  modulus)  of  film  systems  with  the  existence  of  free 
surfaces  at  nanoscale.  Taking  the  calibrated  force-field  parameters  as  a  reference,  we 
systematically  vary  s  from  0.2  to  0.4  kcal/mol  to  generate  polymer  models  with  different 
intermolecular  interaction  strength,  while  keeping  a  the  same  and  thus  the  density  nearly 
constant  for  all  simulations.  It  should  be  noted  that  once  £  is  varied,  it  is  no  longer  representative 
of  the  original  PMMA  CG  model.  The  similar  two-bead  type  of  CG  models  of  methacrylate 
polymers  has  been  recently  reported  Comparing  with  the  two-bead  models,  the  bonded 
interactions  in  our  one-bead  model  preserve  the  general  backbone  features  of  methacrylate 
polymers,  and  the  different  values  of  £  can  be  a  representative  of  different  interactions  resulting 
from  the  different  type  of  side-chain  groups.  The  simpler  single  bead  model  with  less  number  of 
degree  of  freedom  and  atomic  friction  gives  rise  to  overall  softer  mechanical  responses  than 
more  realistic  two-bead  CG  models  and  all-atomistic  models  so  the  study  here  is  meant  to  be 
largely  comparative  and  qualitative. 

Polymer  Free-Standing  Thin  Film  Configuration 

CGMD  free-standing  thin  film  simulations  are  carried  out  using  50  to  500  chains  with  100 
repeat  units  per  chain,  resulting  in  thin  films  with  thicknesses  ranging  from  ~5  to  '-'90  nm  (along 
the  z  axis).  Periodic  boundary  conditions  of  the  simulation  box  are  applied  in  x-y  plane  with 
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dimensions  of  9  nm  x  9  nm  except  for  -'5  nm  thick  films,  in  which  we  use  larger  dimensions  of 
18  nm  X  18  nm  in  x-y  plane  to  improve  the  sampling.  The  film  thickness  is  estimated  as  the 
difference  between  the  maximum  and  minimum  z  coordinate.  In  order  to  create  two  free  surfaces 
of  the  film,  the  box  dimension  along  the  z-axis  are  taken  to  be  larger  than  the  film  thickness  and 
the  boundary  condition  in  that  direction  is  nonperiodic.  The  initial  film  structures  are  built  based 
on  a  self-avoiding  random  walk  algorithm  and  placed  at  the  center  of  the  simulation  box  along 
the  z-dimension.  Then,  the  film  first  undergoes  two  annealing  cycles  from  550  K  to  210  K  and 
then  equilibrated  at  250  K  for  2  ns  under  the  isothermal-isobaric  ensemble  (NPT).  The  time  step 
is  chosen  to  be  4  fs. 

Glass  Transition  Calculation 

The  calculations  of  glass  transition  temperatures  of  both  bulk  and  film  samples  are  carried  out 
using  the  method  described  by  Tsige  and  Taylor  from  the  measurement  of  mean-squared 
displacement  (MSD)  data  of  monomers  in  polymer  chains.  The  MSD  of  polymer  beads  is 
evaluated  as: 

^«W  =  (|fiW-fi(0)|')  (2) 

where  rj(t)  is  the  position  of  the  ith  bead  at  time  t  and  (...)  denotes  the  ensemble  averages.  The 
time-averaged  values  of  Qq  from  40  ps  to  2000  ps  are  plotted  against  temperature  from  490  K  to 
190  K  at  an  interval  of  20  K,  which  provides  a  metric  of  segmental  mobility  as  a  function  of 
temperature.  The  MSD  data  in  the  lower  temperature  region  and  higher  temperature  region  can 
be  fitted  with  linear  slopes,  and  the  intersection  of  the  two  slopes  marks  the  Tg . 


Calculation  of  Elastic  Modulus 
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To  calculate  the  elastie  modulus,  strain-eontrolled  uniaxial  tensile  deformation  is  simulated  in 
the  X  direetion  (along  the  plane  of  the  film),  with  a  eonstant  strain  (engineering  strain)  rate  of  0.5 
ns'\  while  the  pressure  in  the  transverse  direetion  (y)  is  eontrolled  to  be  zero  using  the  Nose- 
Hoover  barostat.  The  z  dimension  of  the  simulation  box  is  automatieally  adjusted  to  shrink  and 
wrap  to  encompass  all  the  beads  during  the  simulations.  The  stress  eomponent  in  tensile 
direetion  is  ealculated  based  on  the  atomie  virial  stress  tensor 


(rai}Urap)j 


V 


a>p  ^ap 


(3) 


where  V  is  the  volume  of  the  actual  polymer  system  instead  of  the  simulation  box,  N  is  the  total 
number  of  CG  beads  or  monomers,  and  and  denote  the  mass  and  veloeity  of  ath 
monomer,  respectively.  The  distance  between  the  monomer  pair  a  and  (3  is  denoted  by  r^p.  U  is 
the  total  potential  energy  including  bond,  angle  and  pair-wise  nonbonded  interaetions.  While  the 
strain  rate  used  in  our  study  is  higher  than  experimental  studies,  it  is  reasonably  low  for  MD 
simulations,  and  the  glassy  modulus  dependenee  on  strain  rate  is  relatively  weaker  than  rubbery 
modulus  measurements  Thus,  qualitatively,  we  don’t  antieipate  strong  rate  dependenee  effeets 
on  the  comparative  analyses  presented  herein.  The  elastic  modulus  is  measured  as  the  slope  from 
0  to  4%  strain  over  10  replieas. 


RESULTS  AND  DISCUSSION 

We  first  characterize  how  Tg  and  elastic  modulus  E  depend  on  the  magnitude  on  the 
intermoleeular  interaetion  strength  (s)  for  bulk  polymer  systems  (no  free  surfaees).  The 
simulation  results  of  bulk  7^  and£'  are  summarized  in  Table  I.  Both  7^  and  £"  systematieally 
inerease  with  £.  Recently,  we  have  shown  that  inereasing  £  from  0.2  (soft)  to  0.4  keal/mol  (hard) 
leads  to  a  monotonie  inerease  in  bulk  7^  whieh  ean  be  aseribed  to  the  faet  that  a  stronger 
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intermolecular  interaction  leads  to  a  greater  activation  energy  barrier  for  polymer  chains  to 
achieve  cooperative  segmental  motion.  The  result  of  uniaxial  tensile  tests  for  glassy  bulk 
polymer  at  T  =  250  K  is  shown  in  Figure  2.  From  the  stress-strain  responses,  there  is  a  distinct 
difference  of  deformation  mechanisms  between  different  £  values.  For  stronger  intermolecular 
interaction  (£  =  0.4  kcal/mol),  four  regimes  can  be  observed:  elastic,  yield,  strain  softening  and 
hardening  regions.  This  is  similar  to  experimental  stress-strain  observations  on  glassy  polymers 
However,  for  weaker  intermolecular  interaction  (f  =  0.2  kcal/mol),  there  is  no  obvious  strain 
softening  or  hardening  regime.  The  strain  softening  and  hardening  regimes  with  a  stronger  £  can 
be  attributed  to  an  increased  activation  energy  associated  with  molecular  arrangement  and 
interchain  sliding  after  elastic  deformation  For  the  elastic  response,  the  bulk  elastic  modulus  E 
increases  with  £,  which  is  consistent  with  previous  findings,  since  the  elastic  work  is  primarily 
stored  as  non-bonded  internal  energy 

In  polymer  thin  films,  regions  near  free  surfaces  characterized  by  an  enhanced  mobility  and  a 
reduced  density  compared  to  interior  bulk-like  regions  due  to  the  less  constraint  on  the  surfaces 
are  often  called  interfacial  layers.  The  thickness  of  the  interfacial  layer,  or  the  so-called 
interfacial  thickness  of  a  free-standing  film  is  calculated  on  the  basis  of  Gibbs  dividing 
surface  (GDS)  concept  in  thermodynamics,  which  can  be  used  as  an  indication  of  free  surface 
effect  on  the  mobility  of  polymer  films.  Figure  3  (a)  shows  the  density  profile  along  the  z 
position  of  the  film.  The  value  of  is  calculated  by  the  equation:  =  (hmax  ~  ^e//)/2, 

where  is  the  maximum  film  thickness  defined  by  the  distance  between  the  maximum  and 
minimum  z  coordinates  of  the  film  and  is  the  effective  thickness  between  the  two  GDS 
positions  as  labeled  by  the  vertical  dotted  lines.  Figure  3  (b)  shows  the  result  of  as  a 
function  of  film  thickness  for  different  £  at  250  K  (below  bulk  Tg).  From  the  result,  it  can  be 
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clearly  observed  that  the  polymer  film  with  a  lower  value  of  s  has  a  larger  and  thus  a  greater 
free  surfaee  effect  on  chain  mobility  for  all  the  film  thieknesses.  The  obtained  interfacial 
thickness  is  about  0.5  to  1.0  nm,  whieh  is  in  the  similar  range  as  reported  by  earlier  MD  studies 
of  polymer  films  For  each  e,  the  film  with  different  thiekness  exhibits  almost  the  same 

value  of  By  taking  the  mean  value  of  over  different  film  thicknesses,  we  obtain  the 
average  interfacial  thickness  for  each  s,  which  is  listed  in  Table  II. 

The  thickness  dependence  of  of  free-standing  films  with  different  intermolecular 

interaetions  is  shown  in  Figure  4.  The  behaviors  that  the  increases  with  film  thickness  and 
then  converges  to  their  bulk  Tg  are  observed  for  all  the  e  values  and  eomparable  with 
experimental  observations.  The  film  exhibits  a  greater  jp^^  over  different  thickness  for  a  larger 
value  of  e.  However,  for  a  polymer  with  a  larger  e,  the  jp^^  converges  to  their  bulk  value  faster 
with  inereasing  thiekness  compared  to  a  softer  polymer  with  a  lower  e.  The  free  surface  induced 
reduetion  in  film  7^,  which  can  be  attributed  to  the  enhanced  mobile  layer  at  a  surfaee  as 
discussed  in  earlier  studies  is  more  significant  for  weaker  intermolecular  interaction. 

Therefore,  the  simulation  result  suggests  that  the  free  surfaee  effect  on  the  depression  of  film  Tg 
is  stronger  for  weaker  intermolecular  interaction. 

To  further  quantify  film  thickness  effeets  on  Tg  scaling  and  its  relation  to  intermoleeular 
interactions  as  shown  in  Figure  4,  we  employ  the  empirical  fitting  formula  introduced  by  Kim  et 
al.  to  describe  Tp^^: 
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where  h  is  the  film  thiekness,  is  obtained  from  the  bulk  simulations  (Table  II),  and  5^  is 
ealled  the  eharaeteristie  thiekness  that  determines  the  saturation  rate  ofTp^^  growth.  The 
dashed  lines  in  the  plot  represent  the  best  fit  of  Eq.  4  to  the  data,  and  the  obtained  result  of  5^  is 
summarized  in  Table  IT  By  employing  this  simple  empirieal  form  that  involves  only  one  fitting 
parameter,  we  are  able  to  more  quantitatively  assess  how  the  relative  free  surfaee  effeets  on  the 
film  ehange  with  different  intermoleeular  forees.  From  Table  II,  it  ean  be  observed  that  the  value 
of  5^  depends  on  the  intermoleeular  interaction  s  and  it  decreases  with  increasing  £,  which  is 
qualitatively  consistent  with  obtained  from  thermodynamic  analysis  of  the  film,  even  though 
the  value  of is  smaller  than  for  cache.  This  result  indicates  that  the  thermodynamic 
analysis  of  free  surface  (GDS)  can  be  well  correlated  with  the  thickness  dependence  of  film  glass 
transition  behavior,  and  suggests  that  both  and  5^relate  to  the  depth  of  the  free  surface  layer 
with  enhanced  mobility  for  the  systems  studied  herein. 

Next,  we  measure  the  elastic  modulus  of  free-standing  films  with  different  £,  as  well  as 
different  thickness  values.  Figure  5  (a)  shows  the  snapshots  of  MD  simulations  during  the  tensile 
tests  of  the  free-standing  film.  The  thickness  dependence  of  the  modulus  of  the  films  at 
250  K  is  shown  in  Figure  5  (b).  Similar  to  the  glass  transition  behaviors  of  the  films,  the  free 
surface  induced  decrease  in  compared  to  the  bulk  value  becomes  more  pronounced  as  the 
thickness  decreases. 

It  appears  from  our  results  that  both  and  are  influenced  by  the  presence  of  the 
free-surface  effects  in  a  free-standing  thin  films.  This  finding  leads  to  the  following  question  that 
we  seek  to  answer  here:  does  the  characteristic  length  scale  at  which  the  deviation  of  E^'-^'^ 
emerges  from  the  bulk  value  correlates  to  that  for  the  glass  transition  temperature? 
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As  a  first  attempt  to  elucidate  this  issue,  here  we  directly  compare  the  free  surface  effects  on 
the  modulus  and  Tg  of  free-standing  thin  films  by  employing  the  functional  form  of  Eq.  4  for 

both  and  Specifically,  we  describe  the  thickness-dependence  of  film  modulus  as: 


7  bulk 


Ef“”'{h)  = 


l  +  ((5^//i) 


(5) 


where  5^  is  the  characteristic  thickness  obtained  from  the  modulus  measurement.  The  best  fit  of 
Eq.  5  is  denoted  by  the  dotted  lines  in  Eigure  5  (b). 

In  essence,  5^  and  5^  can  be  considered  as  the  comparable  figures  of  merit  that  quantify  the 
dependence  of  elastic  modulus  and  glass-transition  temperature  on  film  thickness,  which  is 
directly  related  to  length  scale  where  film  properties  deviate  from  bulk  values.  The  simulated 
values  of  5^  as  a  function  of  intermolecular  interaction  strength  are  listed  in  Table  II.  The  value 
of  5^  increases  from  4  to  10  nm  as  £  decreases  from  0.4  to  0.2  kcal/mol.  This  is  reminiscent  of 
the  observation  that  the  free  surface  effect  on  Tg  is  greater  for  lower  e.  The  values  of  5^  are 
greater  than  those  values  of  5^for  the  same  e,  which  indicates  that  the  convergence  of  the 
modulus  to  the  bulk  value  with  an  increase  in  the  film  thickness  occurs  more  slowly  compared 
with  as  shown  in  Eigure  5  (c). 

Our  observations  over  a  broad  parameter  range  hints  towards  a  general  applicability  of  some  of 
the  concepts  arising  from  previous  experimental  studies.  A  study  by  Torres  et  al.  has  shown 
that  polymethacrylate  films  exhibit  a  much  larger  length  scale  for  observed  changes  in  the 
modulus  at  room  temperature  in  comparison  to  their  apparent  Tg.  Eorrest  et  al.  also  showed 
that  the  high-frequency  mechanical  properties  of  free-standing  PS  film  measured  by  Brillouin 
light  scattering  (BES)  are  not  correlated  with  the  film  Tg  as  film  thickness  is  varied.  Wang  and 
Mckenna  report  that  the  significant  Tg  reduction  for  PS  when  the  film  thickness  is  less  than  20 


15 


nm  using  the  liquid  dewetting  method.  Additionally,  earlier  studies  by  MeKenna’s  group  and 
others  have  shown  that  the  observed  ehanges  in  rheologieal  properties  of  polymer  films 

with  varying  film  thieknesses  eould  not  be  well  eorrelated  with  the  Tg  measurement. 

To  further  interpret  our  results  and  quantify  the  length  seale  of  perturbation  of  free  surfaee,  the 
information  of  distributions  of  the  loeal  stress  and  relaxation  assoeiated  with  segmental  mobility 
of  polymer  ehains  within  the  film  ean  provide  insight  into  the  nature  of  the  surfaee  effeet.  Figure 
6  (a)  shows  the  plot  of  both  spatial  distribution  of  loeal  stress  and  segmental  relaxation  time  as  a 
funetion  of  film  position  z.  The  loeal  stress  t(z)  =  -  +fyy(^)]  is  measured  from  the 

eenter  to  the  free  surfaee  of  the  film  intermediately  after  applying  a  biaxial  step  strain  of  10%  in 
X  and  y  directions  for  £  =  0.2  kcal/mol.  It  can  be  clearly  seen  that  the  unrelaxed  local  stress  is 
zero  at  the  free  surface  boundary,  and  increases  more  or  less  linearly  with  the  depth  into  the  film, 
as  labeled  by  the  solid  slope  in  the  plot.  At  -'7  nm  away  from  the  free  surface,  the  local  stress 
saturates  to  the  value  in  the  interior  region,  which  is  indicated  by  the  horizontal  dashed  line. 
Since  at  a  constant  strain,  the  stress  developed  will  be  related  to  the  local  modulus  of  the  film, 
this  analysis  reveals  that  the  free  surface  effect  that  gives  rise  to  a  linear  gradient  in  the  modulus 
that  near  surfaces  can  penetrate  several  nanometers  into  the  film. 

The  local  relaxation  time  Tq  plotted  in  Figure  6  (a)  can  be  defined  as  the  time  needed  by  the 
mean-squared  displacement  of  CG  beads  within  the  local  region  to  reach  the  bead  size  a  =  6.5  A: 

gg(f  =Tg,Z)  =  (7 

Compared  to  the  stress  distribution,  it  can  be  observed  that  the  length  scale  of  the  relaxation 
gradient  within  the  surface  region  is  slightly  shorter  than  that  of  the  stress  gradient.  Considering 
the  fact  that  the  film  elastic  modulus  converges  to  the  bulk  value  at  a  larger  film  thickness 
compared  to  Tg,  this  finding  may  suggest  that  length  scales  of  film  Tg  and  modulus  do  not  have 
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to  be  exactly  the  same  because  the  length  scales  of  the  gradient  of  the  local  stress  and  mobility  in 
the  surface  region  can  be  different.  It  should  be  noted  that  it  is  very  challenging  to  accurately 
predict  the  spatial  Tg  within  very  small  region  less  than  -'1  nm  using  the  current  MSD  approach. 
This  is  mainly  because  the  segmental  mobility  of  polymer  chains  can  be  very  high  and  the 
monomers  can  diffuse  out  of  the  spatial  region  where  they  initially  stay  into  other  layers  at  a 
temperature  close  or  above  Tg .  However,  the  local  mobility  via  the  measurement  of  local  MSD 
should  be  related  to  Tg,  since  both  properties  are  measured  based  on  the  information  of 
segmental  dynamics  in  our  study. 

On  the  basis  of  the  local  stress  distribution  across  the  film,  we  employ  a  composite  bilayer 
model  as  illustrated  in  Figure  6  (b)  to  explain  the  thickness  dependence  of  The  film  with  a 
total  thickness  of  h  can  be  considered  as  a  composite  with  two  softer  surface  layers  with  a 
thickness  of  and  one  stiffer  bulk-like  interior  region  based  on  the  spatial  stress  distribution. 
The  shape  of  the  profile  in  spatial  stress  is  similar  to  that  in  density,  but  the  values  of  is 
larger  than  the  interfacial  thickness  determined  from  GDS.  Based  on  the  observation  that 
is  nearly  independent  of  film  thickness,  one  can  assume  that  this  linear  gradient  holds  for 
different  film  thickness.  Following  this,  the  surface  layer  has  an  average  modulus  of  = 
E^uik  j2^  As  the  distance  from  the  free  surface  increases  beyond  the  linear  gradient  region,  the 
free  surface  effect  vanishes,  and  the  modulus  can  be  considered  to  be  equal  to  If  the  film 

thickness  h  is  less  than  ,  this  means  that  the  effects  of  the  two  free  surfaces  are  pervasive 

throughout  the  film  and  there  is  no  interior  bulk-like  region.  This  approach  is  qualitatively 
similar  to  the  concept  of  the  glassy  bridge  reinforced  model  that  is  introduced  to  explain  the 
mechanical  reinforcement  of  nanofilled  elastomers.  Then,  becomes  a  linear  function  of  h 
and  is  directly  equal  to  Therefore,  the  can  be  obtained  based  on  the  rule  of  mixtures: 


17 


E 


film 


- ) 

h 

r'bulk 


h  s  2/1™'^ 


h  < 


(7) 


The  prediction  of  Eq.  7  is  represented  by  the  solid  curves  in  Figure  5  (b).  For  thinner  films 
with  thicknesses  below  '-lO  nm,  the  prediction  slightly  underestimates  the  simulation  results, 
which  can  be  attributed  to  the  less  accuracy  in  the  assumption  of  linear  modulus  gradient  in  the 
surface  region  when  the  film  thickness  is  less  or  comparable  to  .  However,  without  any 

input  fitting  parameters,  the  behavior  of  the  thickness  dependence  of  is  well  captured  by 
the  bilayer  composite  model.  The  values  of  determined  from  the  local  stress  distribution 
increase  with  decreasing  e  (Table  II),  which  is  consistent  with  our  previous  empirical  scaling 
analysis  using  Eq.  5.  Even  though  the  values  of  arising  from  the  stress  perturbation  of  free 
surface  are  only  in  the  range  of  4  to  7  nm  for  different  £  studied  herein,  the  resultant  change  in 
modulus  can  be  observed  in  a  film  with  a  thickness  up  to  around  100  nm  that  is  significantly 
larger  compared  to  .  These  values  are  also  comparable  with  experimental  observations.  In 
the  previous  experimental  study  ",  the  thickness  of  the  free-surface  layer  for  PMMA  films  on 
PDMS  is  reported  to  be  about  3.5  nm.  Our  simulation  prediction  of  for  e  =  0.4  kcal/mol  is 
quite  close  to  their  value.  In  a  recent  study  by  Brinson  et  al.  they  show  that  the  interphase  in 
the  small  local  region  near  the  substrate  with  a  higher  local  modulus  results  in  an  altered  elastic 
modulus  of  fdm  up  to  hundreds  of  nanometers,  which  provides  a  support  to  our  observation.  In 
the  simulations  by  Yoshimoto  et  al.  free-standing  polymer  film  exhibits  local  softening 
regions  with  lower  storage  modulus  at  two  free  surfaces  via  the  measurement  of  local  dynamic 
mechanical  properties,  and  the  spatial  modulus  exhibits  an  nearly  linear  gradient  in  the  surface 
region,  which  is  very  similar  to  what  we  observe  herein.  Therefore,  our  finding  physically 
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explains  why  the  approaehes  to  the  bulk  value  as  the  film  thiekness  inereases  by 
employing  the  bilayer  eomposite  model,  and  highlights  the  importanee  in  determining  the  length 
seale  of  interphase  region  arising  from  free  surfaees  in  polymer  films. 


CONCLUSION 

In  this  study,  we  provide  a  eomprehensive  CGMD  study  that  quantifies  how  the  glass-transition 
temperature  and  elastic  modulus  at  glassy  state  of  free-standing  polymer  thin  films  depend  on 
intermolecular  interactions  and  film  thickness.  Both  bulk  and  film  systems  exhibit  a  greater  value 
of  7^  and  modulus  with  similar  backbone  structure  but  stronger  intermolecular  interaction. 
However,  the  free  surface  has  a  stronger  effect  on  Tg  and  modulus  for  the  film  with  a  weaker 
intermolecular  interaction  at  nanoscale.  By  comparing  the  results  of  Tg  and  modulus  of  the  film, 
we  find  that  the  convergence  of  the  modulus  to  the  bulk  value  with  an  increase  in  the  film 
thickness  occurs  more  slowly  compared  with  Tg .  This  longer-range  free  surface  effect  on 
modulus  than  Tg  may  be  further  supported  by  the  observation  of  a  greater  length  scale  of  local 
stress  distribution  in  the  surface  region  than  that  of  local  relaxation  time.  The  observed  decrease 
in  modulus  of  film  exists  up  to  the  thickness  over  -'100  nm.  However,  the  measurement  of  local 
stress  distribution  induced  by  a  step  strain  deformation  reveals  that  there  exists  a  soft  region  with 
a  length  scale  on  the  order  of  several  nanometers,  which  is  much  smaller  than  the  critical  film 
thickness  where  the  film  modulus  converges  to  bulk  value.  The  thickness  dependence  of  film 
modulus  is  well  captured  by  employing  the  concept  of  bilayer  composite  model,  which 
physically  explains  why  the  modulus  increases  with  film  thickness.  Our  finding  establishes  a 
comprehensive  description  of  free-standing  film  elastic  properties,  providing  an  analytical 
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description  that  highlights  the  importance  in  determining  the  length  scale  of  interphase  region 
arising  from  free  surfaces  in  polymer  thin  films. 
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Tables 

TABLE  I.  Comparison  of  glass-transition  temperature  (Tg)  and  elastic  modulus  (E)  of  bulk 
polymer  at  a  glassy  state  (T  =  250  K)  for  different  intermolecular  interaction  strength  (e). 


£  (kcal/mol) 

E  (MPa) 

0.2 

321.5  (±18.7) 

139.2  (±26.3) 

0.3 

388.6  (±10.1) 

288.5  (±35.5) 

0.4 

411.1  (±21.7) 

411.0  (±44.3) 

TABLE  II.  Summary  of  different  resultant  surface  related  properties  that  are  used  to  quantify 
the  free  surface  effects  and  their  relevant  length  scales  as  described  in  the  text  for  different 


intermolecular  interaction  strength  (s). 


£  (kcal/mol) 

(A) 

(nm) 

(nm) 

0.2 

8.6 

3.2 

10.3 

7.0 

0.3 

6.9 

2.1 

6.6 

5.0 

0.4 

5.7 

1.0 

4.2 

4.0 
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Figures  and  Captions 


Free-standins  Film 


Film  thickness  (h)  h  =  W 


Coarse-grained  model 


Free-standing 
polymer  film 


FIG,  1,  (a)  Illustration  of  polymer  nanostructures  with  a  width  in  nanoscale  (W  «  L,  H)  (left) 
and  the  configuration  of  all-atomistic  polymer  free-standing  film  in  simulations  with  periodic 
boundary  conditions  in  x  and  y  axis,  (b)  The  mapping  from  the  all-atom  model  of  poly(methyl 
methacrylate)  (PMMA)  to  the  coarse-grained  model  (left)  and  the  snapshot  of  the  coarse-grained 
free-standing  polymer  film  with  two  free  surfaces  on  the  top  and  bottom  faces  (right). 
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FIG,  2,  Stress-strain  response  for  bulk  polymers  with  different  intermoleeular  interaetions  (f  = 
0.2,  0.3  and  0.4  keal/mol)  under  uniaxial  tension  at  250  K.  The  elastie  modulus  is  determined  by 
fitting  the  data  with  the  slope  from  0  to  4%  strain  over  10  replieates. 
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FIG,  3,  (a)  The  density  profile  along  the  height  of  free-standing  thin  film  (z)  for  e  =  0.3  keal/mol 
at  250  K.  (b)  Interfacial  thickness  of  mobile  free  surface  layer  as  a  function  of  film 
thickness.  The  dashed  lines  indicate  the  mean  value  of  over  different  film  thicknesses. 
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FIG,  4,  The  thickness  dependence  of  the  glass  transition  temperatures  of  free-standing  thin  film 
jfiim  different  a.  The  solid  line  is  calculated  as  the  best  fit  of  Eq.  4,  which  results  in  the  fit 
parameter  of  the  characteristic  thickness  for  . 
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(C) 
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FIG,  5,  (a)  MD  simulation  snapshots  during  the  tensile  deformation  of  the  free-standing  film 
with  a  thickness  of  18  nm.  The  arrows  indicate  the  direction  of  deformation  of  the  film,  (b)  The 
elastic  modulus  of  free-standing  film  as  a  function  of  film  thickness  for  different  e.  The 
dotted  line  is  the  fit  of  empirical  scaling  formula  in  Eq.  5.  The  solid  line  is  prediction  of  Eq.  7 
derived  from  the  bilayer  composite  model,  (c)  Comparison  of  the  empirical  scaling  results  for  the 
normalized  Tg  ^  ^  function  of  film  thickness.  (Inset) 

Comparison  of  normalized  E  for  different  s. 


Surface  layer 

m" 


Interior  bulk 
layer 


Surface  layer  /f 


(b) 

FIG,  6,  (a)  Spatial  distributions  of  the  stress  and  relaxation  time  across  the  film  with  a  thickness 
of  ~78  nm  from  the  center  to  the  free  surface  {e  =  0.2  kcal/mol).  The  local  stress  ~ 
{Txx  +  T^yy)  is  measured  in  a  biaxial  deformation  with  a  step  strain  of  10%  in  x  and  y  directions. 
The  local  relaxation  time  is  determined  based  on  Eq.  6.  The  solid  slope  and  dashed  vertical  lines 
correspond  to  the  linear  gradient  of  the  stress  within  the  surface  layer  and  the  average  stress  and 
relaxation  time  in  the  interior  region  of  the  film,  respectively.  The  thickness  of  the  surface  layer 
j^surf  defined  as  the  thickness  of  the  linear  slope  region  of  the  stress  at  the  free  surface  labeled 


by  the  colored  shadow  in  the  plot,  (b)  The  schematic  of  the  bilayer  composite  model  that  consists 
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of  two  soft  surface  layers  with  a  thiekness  of  and  an  interior  bulk-like  layer  derived  from 
the  loeal  stress  distribution. 
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